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The temperature dependence of the symmetry energy and the symmetry free energy coefBcients 
of atomic nuclei is investigated in a finite temperature Thomas-Fermi framework employing the 
subtraction procedure. A substantial decrement in the symmetry energy coefficient is obtained for 
finite systems, contrary to those seen for infinite nuclear matter at normal and somewhat subnormal 
densities. The effect of the coupling of the surface phonons to the nucleonic motion is also considered; 
this is found to decrease the symmetry energies somewhat at low temperatures. 
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I. INTRODUCTION 



The nuclear symmetry energy reflects the energy cost 
in converting the isospin asymmetric nuclear matter to 
the symmetric one. Conventionally, as part of global 
mass formulas in the framework of the liquid-drop model 
(i.e., ignoring pairing and shell-corrections), the symme- 
try energy Es of an even-even or odd-odd finite nucleus 
with N neutrons and Z protons is given from the cxprcs- 



where 



sionEs = E{N,Z)-E{A/2,A/2) = ^sym-^ 
A = N + Z, E{N, Z) is the energy of the nucleus barring 
the Coulomb part, and Csym is the symmetry energy co- 
efficient. For infinite nuclear systems, the value of esym 
is usually taken in the range of 30-34 MeV [H-Q ■ Apart 
from the evident basic role of symmetry energy in the 
correct description of the nuclear binding along the peri- 
odic table and a broad understanding of the nuclear drip 
lines, it is of seminal importance in giving a proper de- 
scription of the dynamical evolution of the core-collapse 
of a massive star and the associated explosive nucleosyn- 
thesis. A large (small) magnitude of Csym inhibits (accel- 
erates) change of protons to neutrons through electron 
capture [1, 0| • This subtly changes the equation of state 
(EOS) of the dense stellar matter and influences the dy- 
namics of the collapse and explosive phase of a massive 
star. 

Stellar matter in the collapse or bounce phase is warm, 
knowledge about the temperature dependence of the 
symmetry energy is therefore a concern of paramount 
importance. It is also an important ingredient in prop- 
erly accounting for the multiplicity distributions of the 
fragment masses and charges and also their isotopic dis- 
tributions in multifragmentation of hot nuclear systems. 
Temperature dependence of nuclear symmetry energy 
and the symmetry free energy Fs (defined likewise Eg as 



F{N,Z)-FiA/2,A/2) ^ f,y, 



[N-zy 



where F is 



the free energy) for homogeneous infinite nuclear matter 
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has thence been studied; model calculations show that Eg 
for infinite matter decreases with temperature while Fs 
exhibits the opposite temperature dependence Q. The 
temperature dependence becomes more prominent at low 
densities, particularly for Eg. In stellar core-collapse or 
bounce phase, matter is, however, not homogeneous, it 
is nucleated 0, Q and therefore it is important to know 
how the symmetry energies of finite nuclei behave with 
increasing temperature. The present communication is 
aimed to understand this feature. 

Calculations on symmetry energy for finite nuclei in a 
low temperature domain {T < 2 McV )have been done 
earlier by Donati et.al for applications in core-collapse 
supernova simulations. At finite temperature, the sym- 
metry coefficient was found to be somewhat larger than 
for cold nuclei [T =0) hindering electron captures. The 
calculations were schematic, but takes into account cou- 
pling of the nucleons to the dynamical surface phonons 
that results in an increased effective nucleon mass, the 
so called w-mass embodying nonlocality in time [l^l ■ 
This effective mass (m^) decreases with temperature, re- 
sulting in an increase in the symmetry energy. Calcu- 
lations by Dean et.al [ll| in a shell model Monte Carlo 
(SMMC) framework support these findings. These cal- 
culations have also been done in the same limited tem- 
perature range and may suffer from shortcomings in the 
isospin dependence of the schematic residual interaction 
they have used. The symmetry coefficients calculated by 
them are much below the accepted nominal values. 

The present calculations have been done in a broader 
temperature domain (T < 8 MeV) to have a better un- 
derstanding of the temperature dependence of symmetry 
energies in atomic nuclei. Finite temperature Thomas- 
Fermi framework with the subtraction technique has 
been employed for this purpose. Two effective inter- 
actions, namely the modified Seyler-Blanchard (SBM) 
[H, [3 and the SkM* [l^ interactions that give practi- 
cally the same symmetry coefficient (~ 31 MeV) at their 
respective saturation densities have been used in our cal- 
culations. Dynamical changes in the energy mass m^^ 
with temperature is also taken into account following the 
prescription of Refs. [3, [l^l . 
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ical framework is briefly outlined. In Sec. Ill, results 
and discussions are given. The concluding remarks arc 
presented in Sec. IV. 

II. THEORETICAL FRAMEWORK 



TABLE L The parameters of the SBM effective interaction 
(in MeV fm units) 



Ci 


Cu 


a 


h 


d 


AC 


348.5 


829.7 


0.625f 


927.5 


0.879 


1/6 



The methodology employed to calculate the symmetry 
energy and symmetry free energy of finite nuclei as a 
function of temperature is outlined in the following. 



A. Effective interactions 

In describing the hot nucleus in the finite tempera- 
ture Thomas-Fermi (FTTF) approach, we have chosen, 
as already stated, two effective interactions with nearly 

31 McV at the 



sym 



the same symmetry coefficient (e 

respective saturation density for nuclear matter). Both 
these interactions are density dependent and both have 
a quadratic momentum dependence. The SkM* interac- 
tion has common usage and FTTF calculations with this 
interaction [H, [l^, [l^l have been performed before, with 
much success. The Seyler-Blanchard [l^ interaction has 
also been very successfully used by Myers and Swiate- 
cki [iO, [in in the context of the nuclear mass formula. 
The modified version of this interaction (SBM) is found 
to reproduce the ground state binding energies, charge 
rms radii. Giant Monopolc Resonance (GMR) energies 
etc., very well for a host of nuclei from -'^^O to very heavy 
systems [13, [2^ . The properties of pure neutron matter 
calculated with this type of interaction is also seen to be 
in good agreement [23| with those obtained from other 
sophisticated interactions. 

The SBM interaction is given by. 



Veff{r,p,p) = Ci^u vi{r,p) +V2ir,p) , 
^^i = -(l-pVfe')/(ri,r2), 



V2 



p{ri)+p{r2) /(ri,r2) 



with 



/(ri,r2) 



g-|ri-r2|/a 

|ri - r2\/a' 



(1) 



(2) 



The isospin dependence in the interaction is brought 
through the different strength parameters Ci for like- 
pairs (n-n,p-p) and Cu for unlike pairs (n-p). The relative 
separation of the nucleons in configuration and momen- 
tum space are given by r = |ri — and p = |pi — P2|. 
The densities at the sites of the two interacting nucleons 
are given by p(ri) and p(r2). The range of the interac- 
tion is a, b and d are the measures of the momentum and 
density dependence in the interaction and k controls the 
stiffness of the nuclear EOS. 

The five parameters C'l, Cu, a, b, and d for a fixed value 



ing (i) the volume energy coefficient for symmetric nu- 
clear matter {ay = —16.1 MeV), (ii) saturation density 
of normal nuclear matter (taken as po = 0.1533 fm~^, 
corresponding to the radius parameter tq = 1.16 fm), 
(iii) the volume symmetry coefficient esym— 31 MeV, (iv) 
the surface energy coefficient of symmetric nuclear mat- 
ter (os =18.01 MeV), and (v) the energy dependence of 
the real part Vau of the nucleon-nucleus optical potential 
= 0.30). 

The procedure for determining the parameters are 
given in detail in Ref . [l^ . The parameter k is obtained 
by reproducing the giant monopolc resonance energies of 
a large number of nuclei employing the scaling model [2^ . 
The parameters of the interaction are listed in Table. I. 
These parameters differ somewhat from the ones used in 
Ref. |l4l | where the volume symmetry coefficient Cgymwas 
taken as 34 MeV. In the present communication, it is 
taken as 31 MeV, a value commonly used by several au- 
thors with some justification from experimental findings 
[^. The effective /c-niass of the nuclcon rrife coming from 
the momentum dependence of the effective interaction, 
comes out to be 0.62m for symmetric matter, where m is 
the nucleon mass. The isoscalar volume incompressibility 



^■^'oo, symmetry incompressibility Ksym[— (9Po" 



0. 



and the symmetry pressure L[= (3po ^'gp™ )po\ 240, 
— 101 and 59.8 MeV, respectively for the SBM interac- 
tion. The corresponding values for SkM* interaction are 
216.7, -155.9 and 45.8 MeV. 



B. Describing the hot nucleus 

In the FTTF approach, the nucleon density profile for 
the hot nucleus is arrived at self-consistently. The details 
for obtaining the density in this method for the SkM* 
force are given in Ref. [l^; for the SBM interaction, 
they have also been presented earlier [l3,[23- However, 
for completeness, we present the salient features below. 

The total energy of the nucleus is given by 



E 



dr[eK{r) +si{r) +ec{r)] 
= Ek + Ej + Ec 



(3) 



where eii'(r), e7(r) and ec{r) are the kinetic, nuclear in- 
teraction and Coulomb interaction energy densities, re- 
spectively; Ek, E] and Ec are the corresponding total 
energies. The energy densities are given by 



p 



-nr{r,p,T) 



(4) 
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^ / ^-i(|r~r'Up-p'|) 



T 

+V2{\r - r'\,p)}{Cinr{r',p') + C„7i_,(/,p')} 



where 



xririr, p) 



dr'dpdp' 



(5) 



(r) = e^TTpp{r) / drV Pp{r')g{r,r') 



3e 



47r 



(6) 



In these equations, nr{r,p, T) is the nucleon distribution 
function. In the Coulomb energy density, the first and 
second terms are the direct and exchange parts, respec- 
tively. Here, Pp{r) is the proton density profile and 



g(r,r') = 



(r + r') - |r - r'\ 



(7) 



In Eqs. (4) and (5), t is the isospin index; if r refers 
to a proton, — r refers to a neutron and vice versa; the 
distribution function n,- is determined from the min- 
imization of the total thermodynamic potential VL (= 
E -TS -J2r t^rNr) [11. Here S is the total entropy 
of the system that can be obtained from Landau quasi- 
particle approximation as. 



5 = - 



j \nr hinr + (1 — rir) ln(l — tir) drdp, (8) 



and are the nucleon chemical potentials. The distri- 
bution function is given as 



n,(r,p,T) = exp\(^ + + p^V^ + V^^ - pr)] '(9) 

In Eq. (9), V^{r) -\-p'^V^{r) is the single-particle poten- 
tial and V!^ (r) is the rearrangement potential originating 
from the density dependence of the interaction. The full 
expressions for the different components of the single- 
particle potential and are given in Ref. [3l and are 
not repeated here. The effective /c-mass of the nucleon 
TTiT-^kir) can be defined as 



rrir.k 



r 1 

rrir 



2V} 



From Eq. (9), the nucleon density 

/OrCr) = J nr{r,p,T)dp 



(10) 



(11) 



is obtained as, 



A-jr 

A*r{r) = -^[2m.,,(r)T]3/2, 
and Jkiv) is the Fermi integral 



-dx 



(13) 



(14) 



/o 1 + exp{x — 77) 

with the fugacity 77 given as 

rir{r) = [pr-VAr)]/T. (15) 

Here Vr{r) is the effective single-particle (SP) potential 
(Coulomb included), defined as 

Vr{r)^[v?{r) + VHr)+5rMr))/T, (16) 

Vc{r) being the single-particle Coulomb potential, 

Vc{r) = 2eV j dr'r'^pp{r')g{r,r') 

_f!(3.2)l/3 (^)l/3^ (17) 
TT 

From Eqs. (4) and (5), the kinetic and nuclear inter- 
action energy densities can be simplified to 

= ^ E ^i2mrMTf^'J,MVrir)) (18) 



+ T? T.('^rnr,kT)'/'V^\r)Jy,{rjAr))- (19) 



/l3 



For a finite nucleus, the total energy density is the sum of 
Eqs. (6), (18) and (19). For infinite systems, the energy 
densities are space independent. For a particular density 
p, the kinetic and potential energies per nucleon can then 
be defined as ck = £k/p and e/ = £// p, respectively. 

At a large distance from the center, particularly for 
neutrons, Vr is zero and the density becomes p-r ~ 
gMT/T .^^jjjQjj jg 2ero at T = as /Ut is negative. At finite 
T, however, pr is not zero and takes a constant value with 
a nonzero pressure at the surface. This makes the system 
thcrmodynamically unstable. Moreover, the density pro- 
file becomes dependent on the size of the box in which 
calculations are done. This problem is overcome in the 
subtraction procedure [13; HH, where the hot nucleus, as- 
sumed to be a thcrmalized system in equilibrium with a 
surrounding gas simulating the effects of the evaporated 
nucleons, is separated from the embedding environment. 
This method is based on the existence of two solutions 
to the FTTF equations, one corresponding to the liquid 
phase with the surrounding gas (Ig) and the other cor- 
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the hot nucleus in thermodynamic equilibrium is given 
by Pt,i{t) = PT,ig{r) " PT,g(j') and it is independent of 
the choice of the box size. The conservation of nucleon 
number of each species Nr of the hot nucleus gives 



[PrAgir) - pT,gir)]dr ^ Nr. 
The energy E of the nucleus is given by 



En, 



(21) 



where Eig and Eg are the total energies of the liquid- 
gas system and the gas alone. From Eq. (8), the total 
entropy could be recast as [26j . 

where is the single-particle occupation function 

frier, Pr,T) = [1 + CXpiSr - Pr) /T]-\ (23) 

and Qr is the subtracted single-particle level density. It 
is given as [l^] , 



gr{er,T) 



{m'luf'Ner-V[%r) 



4V2 



-(m?_j3/y£,-V|(r)]r2dr, (24) 



where is the single-particle potential, i referring to Ig 
or g. 

The free energy is calculated from F = E — TS. In 
terms of the occupation function, the density in the 
FTTF approximation could be written as 



1 



X f{er,Pr,T)der. (25) 



C. The energy-mass 

The effective mass m* has two components, the fc-mass 
and the energy mass (or the cj-mass); it is defined as 



to( )( ), 

m m 



(26) 



where m is the nucleon mass. 

The temperature-dependent fc-mass is given by the 
momentum-dependent part of the single-particle poten- 
tial as discussed earlier. The w-mass originates from 
the coupling of the single-particle motion with the col- 
lective degrees of freedom. The m^^ varies with posi- 
tion and temperature. A self-consistent calculation of 
the temperature-dependent m^^ is very involved and not 
within the scope of the present work. We have therefore 
taken a phcnomenological form \it\ for such that 



1 - 0.4AV3exJ-r,,^n J^^.(27) 



1 dp{r) 



The temperature T and distance r are measured in MeV- 
fm units, p{0) is the central density of the nucleon dis- 
tribution in the nucleus. The collectivity refers to the 
liquid phase only. The density in the above equation is 
p{r) = pig{r) — Pg{r); A refers to the liquid mass. This 



(20) implies m 



ig 



and 



To avoid the complexities arising from the self- 
consistent calculation of the density profile with the in- 
clusion of w-mass, we have adopted a realistic extension 
of the method given in Ref. [13. This is described in 
some detail in Ref. [2(| . 

The subtracted level density corresponding to Eq.(24) 
then modifies to 



9t(£t,T) 



4^/2 



ra 



-{m^r,k)'^'^f^^^^ydr, 



(28) 



the densities in the Ig or g phase modify accordingly, 

l'',l3/2 



1 



27r2;i3 



2m 



T,k' 



Sr-Vl- 



X f{Sr,Pr,T)der- 



(29) 



The chemical potential pr has modified to pr to conserve 
the particle numbers in the nucleus 



Nr = 



J 9t {£t , T)fr {Sr , Pr , T)der ■ 



(30) 



D. The symmetry coefficients 

For nuclear matter, the symmetry coefficient at density 
p can be defined as, 

esy.n{T) = [e(p, X, T) - e(p, X = 0, T)]/X\ (31) 

where e is the energy per nucleon and X is the asymmetry 
parameter X = (p„ — pp)/p). Here p„ and pp arc the 
neutron and proton densities (p = p„ + Pp).Thc kinetic 
and potential parts of the symmetry coefficients e,,y,„ (k) 
and esym(v) can likewise be defined replacing e in the 
square brackets in Eq. (31) by ck and ej, respectively. 

For a nucleus of mass A, the symmetry coefficient, in 
the spirit of the liquid-drop model, can also be defined 
as 

esym{T) = [e(7V, Z, T) - e{A/2, A/2, T)]/X\ (32) 

Here the asymmetry parameter \sX = (N—Z) /A. In this 
equation, e is the energy per particle of the nucleus bar- 
ring the Coulomb contribution. For relatively heavy nu- 
clei, stable systems are usually isospin-asymmetric, there 
the definition given by Eq.(32) may not be operative. In 
such systems, esym could be calculated from 



r _ / /i TT- m\ 
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where Xi and X2 are the asymmetry parameters of the 
nuclear pair. Similar definitions follow for symmetry free 
energy coefficient fsym on replacing energy e by the free 
energy per particle /. 

The value of the symmetry coefficient obtained from 
Eqs. (32) and (33) depends on the choice of the nuclear 
pair, its value is therefore not unambiguous for a partic- 
ular nucleus. However, the local density approximation 
(LDA) could be used to define the symmetry coefficient 
for a specific nucleus once its density profile is known. In 
the LDA, it is given as [l^] 



N - Z\2 



A 



Pn{r) - Pp{r)^^ 



p{r) 



dr. (34) 



Here e^^[p(r),T] is the symmetry energy coefScient of 
infinite matter at temperature T at a value of the density 
p{r). The local isospin density is given by [pn(r) — Pp(r)]. 



III. RESULTS AND DISCUSSIONS 

As already stated, for the microscopic finite temper- 
ature Thomas-Fermi calculations, we have chosen the 
SkM* and the SBM effective interactions. Both these 
interactions reproduce the ground state properties of 
atomic nuclei and also of nuclear matter fairly well. For 
normal nuclear matter at T=0, these interactions have 
symmetry coefhcient esym ^ 31.0 MeV, a commonly ac- 
cepted value with some support from experimental anal- 
yses . Below the saturation density, both these interac- 
tions display nearly the same density dependence for esym 
as shown in Fig. 1. The symmetry coefficient obtained 
from the SBM interaction (given by the red line) looks 
a little stiffer beyond the saturation density. It would 
be interesting to compare the symmetry coefficients of 
cold nuclei obtained from these interactions with phe- 
nomenological values. This is shown in panels (a) and 
(b) of Fig. 2; the full green lines in both the panels cor- 
respond to the upper and lower bounds for Csym as given 
in Ref. [l as 



1 



04-1/3 ' 



(35) 



where 31.0 MeV and | =2.4 ±0.4. The black 
circles and red triangles refer to calculations with 
SkM* and SBM interactions, respectively. The upper 
panel corresponds to calculations done in the difference 
method using Eq. (33). The nuclear pairs {A, Zi) 
and {A, Z2) collectively written as (A, ^1,^2) are 
chosen as (26,10,12), (40,16,18), (56,24,26), (80,34,36), 
(120,50,52), (150,60,62), (197,77,79),and (238,90,92). 
These results are representative, the values of esym 
depend on the choice of nuclear pairs as shown later. 
The lower panel displays results calculated in the LDA 
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FIG. 1: (color online)The symmetry energy coefficient esym 
for nuclear matter at T =0 as a function of density with SkM* 
(black) and SBM (red) interactions. 



the /3-stability line; they are (26,12), (40,18), (56,26), 
(80,34), (120,50), (150,62), (197,79), and (238,92), 
respectively. In the LDA, the values of Cgym are seen 
to have a relatively weak dependence on the choice 
of the asymmetry parameter X. Compared to those 
in SBM, the SkM* values are somewhat lower. The 
values of Cgym., however, are seen to lie within the 
phenomenological limits, particularly in the LDA. The 
difference in results in Csym for finite nuclei in spite of 
the apparent similarity in its density dependence for 
infinite nuclear matter stems from the fact that the 
SBM is a finite range interaction as opposed to SkM* 
(zero range) which results in somewhat different density 
profiles of nuclei. 

Before embarking on the evaluation of the tempera- 
ture dependence of the symmetry coefficients for finite 
systems, it may be worthwhile to see how Csym or fsym 
behave with increasing temperature for infinite nuclear 
matter. At saturation density, we find that both Csym 
and fsym , in the temperature range we have studied, have 
a very weak temperature dependence with both the inter- 
actions chosen. At relatively low densities, esym is seen 
to decrease with temperature, in consonance with that 
reported earlier by Xu et.al., @ and Moustakidis [28j . 
At these densities, fsym, however, displays a compara- 
bly more prominent rise, in fair agreement with those 
obtained earlier 0, This is displayed in Fig. 3 at 

p = po/8 {po = saturation density) for the SBM inter- 
action. Results with SkM* interaction do not show any 
different behavior and arc therefore not shown. The fall 
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FIG. 2: (color online)The dependence of symmetry energy 
coefficient esym on nuclear mass number with SkM* (black 
circles) and SBM (red triangles) interactions. The limits in 
the nuclear mass formula as given in Ref. [2| are also shown 
(green lines). The upper and lower panels correspond to cal- 
culations using Eqs. (33) and (34), respectively. For details, 
see text. 
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FIG. 3: (color online) In the upper panel the temperature 
evolution of kinetic (green) and potential (black) components 
esym{k) and esym{v), respectively, of symmetry energy coef- 
ficient along with their sum Ssym (red) are shown for nuclear 
matter at a density po/8 calculated with SBM interaction. In 
the lower panel, the same is shown for the total symmetry 
free energy coefficient (red) fsym- 



low density is essentially due to the decrement in its ki- 
netic energy part esym{k) (shown by the green line). This 
part of the symmetry energy decreases at lower densities 
as Pauli-blocking becomes less important because of in- 
creased difFuscncss of nucleon Fermi surfaces with tem- 
perature. The potential component esym{v) (black line) 
decreases more slowly. The increase of fsym with tem- 
perature can be understood from the fact that 



-sym 



(36) 



where s.s„m, the symmetry entropy is negative and is 
~ [23 for dilute nuclear matter. Thus even if 

^sym falls slower, the increase in fsym (lower panel) is 
noticeable. 

In Fig. 4, the temperature dependence of the symme- 
try energy coefficients and their kinetic and potential 
components for a nucleus ^=56 is shown for the SBM 
interaction. The upper panel represents calculations in 
the difference method where Eq. (33) is exploited to ob- 
tain the symmetry coefficients with Zi =24 and Z2 =26. 
The density profiles of the nuclei have been calculated 
with Coulomb switched on, however, the symmetry co- 
efficients are calculated from the energies obtained af- 
ter subtraction of the Coulomb contribution. The lower 
panel shows the calculations in LDA with the use of 
Eq. (34) for the nuclei with Z =24 and 26. The difference 
method shows that for nuclei, the main contribution to 
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FIG. 4: (color online) The upper panel represents the temper- 
ature evolution of kinetic (green) and potential (black) com- 
ponents of symmetry energy coefRcient along with their sum 
esym (red) obtained from A=56 isobaric pair with Z=24 and 
26 using SBM interaction. In the lower panel, the same is 
shown for these isobaric nuclei individually in the local den- 
sity approximation. 
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FIG. 5: (color online) The temperature dependence of 
Esym using SkM* interaction with (SkM*(l)) and without 
(SkM*(0)) inclusion of w-mass are shown in panel (a) by black 
full and dash lines, respectively, for A=56 calculated with 
Z=(24,26) isobaric pair. The red lines represent the same for 
SBM interaction. The same is shown in panel (b) for A=112 
using Z=(48,50) isobaric pair. The panels (c) and (d) repre- 
sent the same for fsym. 
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FIG. 6: (color online) The temperature dependence of esym 
(upper panel) and fsym (lower panel) for the nucleus A=56. 
The black and red lines correspond to results with SkM* and 
SBM interaction, respectively. The numbers in the parenthe- 
ses refer to charge numbers of the different isobaric pairs used 
to calculate the symmetry energies. 



FIG. 7: (color online) The same as in Fig. 6 for the nucleus 
A=64 with isobar pairs as shown in the parentheses. 



part (the black line), esym{k) (the green line) is quite 
small. The symmetry coefScient is nearly constant upto 
~ 2 MeV, but then falls substantially with temperature. 
The calculations in the LDA, for the isobars with Z = 24 
and 26, however, show a very weak temperature depen- 
dence for the symmetry coefficients, which comes mostly 
from the surface region; the temperature dependence of 
infinite matter is refiected in these calculations for finite 
nuclei. Here the contributions from esym{v) and Csymik) 
are comparable. 

Till now, calculations have been reported without the 
inclusion of the effects of the w-mass. As can be seen 
from Eq. (27), nii^/m may be significantly larger than 
unity, more so for heavier nuclei and at lower tempera- 
tures. This would decrease the effective kinetic energy 
of the system and may also decrease the corresponding 
symmetry energy. The effects of w-mass on the symme- 
try coefficients Csym said fsym sltc shown for both the 
SBM (red lines) and SkM* (black lines) interactions in 
Fig. 5. The number (1) in parentheses refers to calcula- 
tions (full lines) with inclusion of the effects of w-mass, 
(0) refers to calculations (dashed lines) without it. The 
results from the difference method are only reported. A 
relatively light {A =56) and a medium-heavy {A =112) 
nucleus are chosen for this purpose. As seen from the up- 
per and lower panels in the figure, Csym and fsym decrease 
with the inclusion of w-mass effects. The dilution of the 
effects due to energy mass with increase in temperature 
is evident from the merging of the full and dashed lines 
with the same color representing a given interaction at 
high T. As stated earlier, the reduction in the symmetry 
coefficients with inclusion of oj-mass is more prominent 
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FIG. 8: (color online) The same as in Fig. 6 for the nucleus 
A=112. 



The dependence of the nuclear symmetry energies on 
the choice of nuclear pairs for a given isobar is shown in 
Fig. 6 for A = 56. The chosen nuclear pairs {Zi, Zi) are 
(24,26), (24,28), and (26,28). The upper panel displays 
the temperature dependence of esy,„ and the lower panel 
shows the same for ] sym ■ The figure shows the sensitivity 
of the symmetry coefficients on the choice of the nuclear 
pair; however, for all the cases, esym decreases whereas 
jsyra increases with temperature. With the nuclear pair 
Z\ = 26 and Z2 = 28, esym even becomes negative at 
high temperature. This violates the general understand- 
ing of the symmetry energy; this negativity is likely to 
arise from the Coulomb polarization of the density pro- 
files of nuclei. We have checked that calculations with 
Coulomb switched off never give negative symmetry en- 
ergy coefficients. The different interactions SkM* and 
SBM yield somewhat different values for egym and fsym', 
in general, values obtained for the symmetry coefficients 
with the SkM* interaction are lower compared to those 
for SBM interaction at lower temperatures. A crossover, 
however, is generally seen at higher temperatures. Cal- 
culations have been performed for nuclear isobars with 
A = 64: and A = 112 also. They are reported in Fig. 7 
and Fig. 8. These results follow nearly the same trend as 
those for A = 56. 



IV. CONCLUDING REMARKS 

We have performed calculations for the temperature 
dependence of the symmetry energy and symmetry free 
energy coefficients of atomic nuclei in a Thomas-Fermi 
model. Calculations have been done with two effective in- 
teractions, namely, the SkM* and the SBM. The subtrac- 
tion procedure has been employed for modeling the hot 
nucleus. Effects of the coupling of the surface phonons to 
the intrinsic particle motion have been taken into account 
through a phenomenological parametric form of the nu- 
cleon cj-mass. The w-mass increases the nucleon effective 
mass at relatively low temperatures, this has an appre- 
ciable effect on the symmetry coefficients which decrease 
noticeably once effects due to nii^ are taken into account. 

Questions may arise about the proper definition of the 
symmetry coefficients for finite nuclei. We have taken 
two definitions, one derived from the usual difference 
method [llj, another from the local density approxima- 
tion (LDA). In both approaches, we find that egym de- 
creases with temperature, fsym shows the opposite trend. 
In the LDA, the fall in Cgym is slow, in the difference 
method, this fall with temperature is much stronger. In 
this method, the results, however, depend sensitively on 
the choice of the nuclear pairs taken to evaluate the sym- 
metry energies. Occasionally, the symmetry coefficients 
may even be negative, this is violative of the general con- 
cept of the symmetry coefficients. 

The ambiguities arising from the different definitions 
for symmetry coefficients of finite nuclear systems may 
possibly be overcome by calculating the ground state en- 
ergies of a host of nuclei over the periodic table with 
a suitable choice of interaction in a microscopic frame- 
work, such as Hartree-Fock (HF) and then finding the 
mass parameters a„ , ag etc. in the spirit of the liquid- 
drop model. The calculations can be extended to various 
temperatures employing finite temperature HF method 
and then find the temperature dependence of the mass 
parameters which is likely to give ambiguity-free tem- 
perature dependence of the symmetry coefficients. This 
is, however, an ambitious program and the work is in 
progress. 
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